Astronomy & Astrophysics manuscript no. 09454yu 


© ESO 2009 


August 24, 2009 





Hot subdwarfs from the stable Roche lobe overflow channel 

S. , and L. Li^ 



o 
o 

(N 

^' 

O 



(N 
> 

m 
o 

On 
O 



' Armagh Observatory, College Hill, Armagh BT61 9DG, Northern Ireland, U.K. 
e-mail: syuOarm . ac . uk 

- National Astronomical Observatories / Yunnan Observatory, the Chinese Academy of Sciences, P.O.Box 110, Kunming, 650011, 
China 



Preprint online version: August 24, 2009 



ABSTRACT 



Context. Hot subdwarfs are core-helium-buming stars with extremely thin envelopes. We discuss the formation and evolution of hot 
subdwarfs formed through the stable Roche lobe overflow (RLOF) channel of intermediate-mass binaries, although their formation 
channels are various. 

Aims. In this study, we concentrate on the formation and evolution of hot subdwarfs binaries through the stable RLOF channel of 
intermediate-mass binaries. We aim at setting out the properties of hot subdwarfs and their progenitors, so that we can understand the 
formation and evolution of hot subdwarfs. 

Methods. Employing Eggleton's stellar evolution code, we have computed conservative and nonconservative population I binary 
evolution sequences. The initial mass of the primary ranges from 2.2 to 6.3 Mq, spaced by approximately 0.1 in logM, the initial 
mass ratio q\ = MJMi is between I.I and 4.5, and the Roche lobe overflow begins at the main sequence, the Hertzsprung gap and 
the first giant branch. In nonconservative binary evolution, we assume that 50 percent of the mass lost from the primary leaves the 
system, carrying away the specific angular momentum of the primary, and the remaining mass is accreted on to the secondary during 
the RLOF. Also, we have studied the distributions of the mass and orbital periods of hot subdwarfs using the population synthesis 
approach. 

Results. We have obtained the ranges of the initial parameters of progenitor binaries and the properties of hot subdwarfs through 
the stable RLOF channel of intermediate-mass binaries, e.g. mass, envelope mass and age of hot subdwarfs. We have found that hot 
subdwarfs could be formed through stable Roche lobe overflow at the main sequence and Hertzsprung gap. We have also found that 

some subdwarf B or OB stars have anomalously high mass (~1 Mq) with a thick envelope (~0.07 Mg 0.16) in our models. By 

comparing our theoretical results with observations on the hot subdwarfs in open clusters, we suggest that more hot subdwarfs in 
binary systems might be found in open clusters in the future. 

Key words, stars: subdwarfs - stars: horizontal branch - stars: evolution - stars: formation. 



1. Introduction 

Hot subdwarfs are generally considered to be core-helium- 
burning stars with extremely thin hydrogen-rich or helium- 
rich envelopes. Earlier observations (Humason & Zwicky [19471 
Feige [T9581 Greens tein [T9661 Slettebak et al. [iMTl Klemola 
11962 ', Berger '1963 ) showed that a number of faint blue field 
stars with early type at high galactic latitudes display peculiar 
spectra. Sargent & Searle (1968 ) defined subdwarf B (sdB) stars 
as stars with colors corresponding to those of B stars in which 
the Balmer lines are abnormally broad compared to population I 
main sequence B stars and they have weaker Hel lines for their 
color Similarly, they defined subdwarf O (sdO) stars as stars that 
have strong Balmer lines relative to main sequence stars and in 
which Hell /14686 is seen. An intermediate class (subdwarf OB, 
sdOB stars), which has effective temperature and photospheric 
helium abundances between those of sdB stars and sdO stars, 
has been reported by Baschek & Norris Cl975i and Hunger et al. 
( I1981I I. As a blueward extension of the horizontal branch (HB), 
hot subdwarfs include three subgroups from their spectroscopic 
classification (Heber 1986): 
(i) sdO - display strong Hell or Hel lines. 
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(ii) sdOB - display hydrogen lines and helium lines. 

(iii) sdB - display strong broadended hydrogen lines and weak 

helium lines. 

The origin of hot subdwarfs remains unclear. In the past two 
decades, two scenarios for the origin of hot subdwarfs have been 
proposed, i.e. a single star model and a binary model. The single 
star model was prevalent in the 1990s because of the absence of 
confirmed hot subdwarfs binaries. The model assumes that the 
envelope of a giant may be striped off by stellar wind near the 
tip of the first giant branch (FGB) and leave an almost bare he- 
lium core. If the helium core is ignited, the star may become a 
single hot subdwarf (Lee 1 19941 Yi, Demarque & Oemler [19971 
D'Cruz et al. 1996). Sweigart (1997) has studied the evolution 
of globular-cluster stars and suggested that helium mixing driven 
by internal rotation substantially increases the helium abundance 
in the envelope. This may result in enhanced mass loss along the 
FGB and the formation of a hot subdwarf. However, the single 
star model is inconsistent with the latest observations of hot sub- 
dwarfs and the UV upturn of elliptical galaxies which is a strange 
phenomenon of the rising flux with decreasing wavelength be- 
tween the Lyman limit and 2500 Afound in almost all elliptical 
galaxies. (Yi & Yoon 2004; Han, Podsiadlowski & Lynas-Gray 
120071) . Recent observations indicate that a quantity of hot sub- 
dwarfs are found in binaries. Maxted et al. ( |2001| l argued that 
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more than two-thirds of their candidates are binaries with short 
orbital periods from hours to days, and Napiwotzki et al. ( |2004| l 
have found that roughly two-fifths of their sample are hot subd- 
warfs binaries. In order to explain the formation of hot subdwarfs 
in binary, Han et al. (I2002I I suggested a binary model. In such a 
model, three main channels are responsible for the formation of 
hot subdwarfs i.e. the common-envelope (CE) ejection channel 
for hot subdwarf binaries with short orbital periods (0.05 ~ 40 
days), the stable Roche lobe overflow (RLOF) channel for hot 
subdwarf binaries with long orbital periods (0.5 ~ 2000 days), 
and the double helium white dwarf (WD) merger channel for sin- 
gle hot subdwarfs. In the CE ejection channel, the progenitor of 
a hot subdwarf experiences dynamical mass transfer on the FGB 
which would result in a CE and a spiral-in phase, leaving a very 
close binary after the envelope has been ejected. If the helium 
core of the giant was ignited, it would become a hot subdwarf 
in a short-period binary with a companion of a main sequence 
or white dwarf. In the stable RLOF channel, the progenitor of 
the hot subdwarf would also lose most of its envelope to pro- 
duce a hot subdwarf through mass transfer, but the mass transfer 
would be stable in contrast to the CE ejection channel. In the 
merger channel, two helium white dwarfs in close binary would 
be driven together by the orbital angular momentum loss due to 
gravitational wave radiation or magnetic braking effects. When 
the white dwarfs merge and the merged object ignites helium, 
this would produce a single hot subdwarf (Webbink 1984; Iben 
& Tutukov |I986l Han 1998). The binary model has explained 
the observations and the UV upturn of elliptical galaxies (Han, 
Podsiadlowski & Lvnas-Grav 120071) . 

The majority of field hot subdwarfs found by recent ob- 
servations (Maxt ed et al. JOOl, Morales-Rueda et al. 120031 
Napi wotzk i et al. 120041 Edelmann et al. [2005 ; Morales-Rueda 
et al. 120051 1 is in short-period binaries, while few hot subdwarfs 
in close binaries were found in globular clusters (Moni Bidin et 
al. '2006a, 2006b). This implies that there would be different for- 
mation channels for field and cluster hot subdwarfs. Burleigh 
et al. d 19991 1 suggested a blue star, which was discovered by 
Elson et al. ('1998) in the young cluster NGC 1818 in the Large 
Magellanic Cloud (LMC), is a possible hot subdwarf. If hot sub- 
dwarfs could exist in very young clusters, their formation would 
be a puzzle. 

Various channels for the formation of hot subdwarfs have 
been studied in the binary model (Han et al. 2002), as well as the 
stable RLOF channel by which a low-mass giant loses most of its 
envelope on the FGB. The giant would leave a degenerate core 
after the mass transfer stops. If the mass of the degenerate core 
is high enough, it would experience a helium flash and the star 
would appear as a core-helium-burning hot subdwarf in a binary. 
In their model, binary evolution sequences with the donor mass 
in the range from 0.8 to 1.9 Mq with the accreting components 
of WD stars or neutron stars (NS) were calculated. In addition, 
it was assumed that mass transfer takes place at the FGB, and 
all the mass lost from the system carried away its orbital angu- 
lar momentum of the accreting component. For the purpose of 
perfecting the stable RLOF channel and attaining a comprehen- 
sive understanding of the formation of hot subdwarfs following 
Han's calculation, we have executed a computation of evolution 
sequences of intermediate-mass binaries. 

In this study, we have simulated the formation and evolution 
of hot subdwarfs through the stable RLOF channel in the con- 
servative and nonconservative case. We have taken into account 
the intermediate mass star filling its Roche lobe at the main se- 
quence (MS), the Hertzsprung gap (HG) or the first giant branch 
(FGB). The mass transfer stops once the radius of the donor is 



smaller than its Roche lobe. If the mass of the helium core of the 
donor (the primary. Mi) is high enough, helium is ignited and the 
star becomes a hot subdwarf in a binary in which the other com- 
ponent is probably a main sequence star or a subgiant star. We 
have obtained the properties of the hot subdwarfs and the initial 
parameters space of their progenitors, e.g. the effective temper- 
ature, the surface gravity, mass, envelope mass, lifetime of hot 
subdwarfs phase. Subsequently, we carried out a Monte Carlo 
simulation in order to acquire distributions of properties of hot 
sudwarfs from the stable RLOF channel of intermediate-mass 
binaries and compared these distributions with observations. 

The outline of this paper is as follows. In section 2, we de- 
scribe the stellar evolution code adopted in this study. In Section 
3, we present the results. The results are discussed in Section 4 
and summarized in Section 5. 

2. The stellar evolution code and the binary 
population synthesis code 

2.1. The stellar evolution code and stellar models 

We use Eggleton's (119711 119721 11973I I stellar evolution code, 
which has been updated with the latest physics over the last three 
decades (Han, Podsiadlowski & Eggleton [T994l Pols et al. 1 19951 
1998). The code uses a self-adaptive non-Lagrangian mesh and 
both convective and semiconvective mixing are treated as a diffu- 
sion process. The stellar structure equations, the mesh equation 
and the chemical composition equations are solved simultane- 
ously. 

The current code uses an equation of state that includes pres- 
sure ionization and Coulomb interaction (Pols et al. 1995) , the 
latest opacity tables derived from Iglesias & Rogers d 19961 1 and 
Alexander & Ferguson (|1994ai |1994bl l via quadratic interpola- 
tion for X = 0.8, 0.7, 0.5, 0.35, 0.2, 0.1, and for F = 0.5 - 
Z, 0.2 - Z and 0. Nuclear reaction rates come from Caughlan 
& Fowler ( 119881 ) and Caughlan et al. (,1985i) , and neutrino loss 
rates are from Itoh et al. ( fT989l[T992b . 

We use a typical Population I composition with hydrogen 
abundance X = 0.700, helium abundance Y - 0.280 and metal- 
licity Z = 0.020 in our computations. We set a = l/Hp, the ratio 
of the mixing length to the local pressure scaleheight is equal 
to 2. Such a value of a gives a roughly correct lower main se- 
quence, as determined observationally by Andersen ( |1991| l. It 
also well reproduces the location of the red giant branch in the 
Hertzsprung-Russell (HR) diagram for stars in the Hyades su- 
percluster (Eggen il985l) , as determined by Bessell et al. (I1989I I. 
A fit to the Sun also result in a = 2 as the most appropriate choice 
(Pols etaL [T998t . 

Convective overshooting is important for the remnant mass 
of a binary after RLOF evolution because the overshooting di- 
rectly affects the scale of the nuclear reaction region in a star. 
In this paper, we follow the work of Schroder, Pols & Eggleton 
( 119971) and use an approach based on the stability criterion itself, 
the 6ov prescription, by incorporating a condition that mixing oc- 
curs in a region with V,- > Va - 5, with 6 defined as the product of 
a specified constant 6ov, the overshooting parameter and a con- 
veniently chosen factor that depends only on the ratio ^ of the 
radiation pressure to the gas pressure: 

6 = 5ov/(2.5 + 20^ + 16^2). (1) 

We take 6ov = 0.12, which leads to overshooting lengths /ov be- 
tween 0.25 and 0.32//p for the mass range of 2.5~7 Mq with 
typical 6 values from 0.05 to 0.04. Such dov best fits eclipsing 
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binaries (Schroder, Pols & Eggleton [T997l Pols et al. [T997] l. and 
here we assume that the influence of convective overshooting is 
the same in our binaries as in the eclipsing binaries. 

RLOF and stellar wind are involved in our models. RLOF is 
treated as a modification of a surface boundary condition, which 
is written as: 



dm ntar T. 

— -CMax[0,(— )^], 
at riobe 



(2) 



where dm/dt is the mass changing rate of the star, rjtar is the 
radius of the star and riobe is the radius of its Roche lobe. In prin- 
ciple, the radius of a primary would equal the radius of its Roche 
lobe if the C is large enough. If the C is very large, however, the 
Roche lobe overflow would not be continuous. Here we take C - 
lOOOMoyr ' following Han et al. ( |2002] l. so that RLOF proceeds 
steadily and the lobe-filling star overfills its Roche lobe as nec- 
essary but never overfills its lobe by much [typically {r^uT /f\obe 
- 1) < 0.001] (Han et al. 1995"2000"2002 ). For steUar wind, we 
adopt Reimers' wind (Reimers .1975J mass-loss law: 



Mwind - 4 X 10-"7//?L/M, 



(3) 



where we use a coeflicient rj - 1/4 (Carraro et al. 119961 Iben & 
Renzini 1983 Renzini 1981). 

Our models include two cases for binary evolution: 

(a) conservative case; 

(b) nonconservative case. 

We assume -AM] is lost from the primary, -/3aMi is accreted 
on to the secondary, then -(1 -yS)AMi is lost from the system, 
carrying away the same specific angular momentum as the center 
of mass of the primary. The change of the angular momentum of 
the system Aj is 



A/ _ (1 -/3)AMiM2 

~r ~ Ml (Ml -I- ' 



(4) 



where J is the orbital angular momentum of the system. Mi is 
the mass of the primary, and M2 is the mass of the secondary. 
We take the mass transfer efficiency fi - l.Q for the conservative 
case, and /3 = 0.5 for the nonconservative case, since the results 
of Paczynski & Zidlkowski ( 1967) and Refsdal, Roth & Weigert 
( 119741 ) indicate that /3 is around 0.5 in the nonconservative case. 
The parameter space for the model grid is three dimensional: 

(i) The range of initial mass of the primaries (Mn): 

logMii = 0.35 (2.23M0) ~ 0.80 (6.3 IMq), and AlogMn = 
0.05; 

(ii) The initial mass ratio q-, - Mi /M2: 
qi = 1.1, 1.5,2.0,3.0,4.0,4.2,4.5; 

(iii) The initial orbital periods of the binaries from the minimum 
period, at which a zero-age main sequence (ZAMS) star 
would fill its Roche lobe, to the maximum period, at which 
a star would fill its Roche lobe at the FGB. 

In this paper, we focus just on the stable RLOF chan- 
nel and do not discuss the case of dynamically unstable mass 
transfer. In each calculation, we checked whether mass transfer 
was dynamically stable. If the mass transfer rate keeps increas- 
ing rapidly, dynamical instability would occur as discussed by 
Podsiadlowski et al. (2002). This leads to the stellar evolution 
code breaking down. If mass transfer is stable, we continue mass 
transfer until the donor shrinks below its Roche lobe, terminating 
the mass transfer phase. In the case of stable mass transfer, if the 
mass of the helium core exceeds the appropriate minimum core 
mass for helium ignition, it will appear as a core-helium-burning 
hot subdwarf. 



2.2. Monte Carlo simulation parameters 

In order to obtain distributions of properties of hot subdwarfs 
from the stable RLOF channel of intermediate-mass binaries, we 
have performed a simple Monte Carlo simulation where we fol- 
low the evolution of a sample of a million binaries. The physical 
inputs of the simulation are depicted as follows: 

(i) We assume that the star formation rate (SFR) is constant 
over the last 13.5 Gyr. 

(ii) For the initial mass function (IMF) of the primary, we 
adopted the IMF of Miller & Scalo (119791 ); the primary mass is 
generated using the formula of Eggleton, Fitchett & Tout 819891 ). 



Ml 



0.19X 



(1 _X)0.75 + 0.032(1 



(5) 



where X is a random number uniformly distributed between 
and 1. The studies by Kroupa, Tout & Gilmore (119931 ) and 
Zoccali et al. (2000) support this IMF. 

(iii) For the initial mass ratio distribution, we suppose a con- 
stant one (Mazeh et al. |T992| Goldberg & MazehUS!!]), 



nil/q) = 1,0 < l/q < 1, 



(6) 



where qi = M1/M2. 

(iv) For the distribution of initial orbital separations, we em- 
ployed that used by Han et al. (2003), where they assume that 
all stars are members of binary systems and that the distribution 
of separations is constant in loga (a is the separation) for wide 
binaries and falls off smoothly at close separations: 



an{a) = 



^ asep,ao < a< ai. 



(V) 



where asep~0.070, oq = IORq, oi = 5.75 x IQ^Rq = 0.13 pc, 
m X 1.2. This distribution implies that there is an equal number 
of wide binary systems per logarithmic interval, and that approx- 
imately 50 per cent of stellar systems are binary systems with 
orbital periods less than 100 yr. 

The total numbers of observed hot subdwarfs are associated 
with their evolution lifetime. Each binary system in our models 
is multiplied by a factor of jt, where yt - t/2 x 10*^ yrs, t is the 
lifetime of hot subdwarf phase. 

3. Results 

In general, three subgroups are defined from photometric and 
spectroscopic classifications: sdB, sdOB and sdO. Fig.l displays 
the position of observed hot subdwarfs. It is seen in Fig.l that 
there is a boundary between sdB, sdOB and sdO stars at effec- 
tive temperature Teif «i35,000 K. The effective temperature of 
sdB and sdOB stars is between 20,000 K and 35,000 K while 
that of sdO is between 35,000 K and 80,000 K. Due to the tight 
correlation of spectral type (effective temperature, T^ff) and spec- 
troscopy, we define three subgroups from spectral type {Tgff) 
classification corresponding to the spectroscopic classification: 

(1) sdB - Teff is between 20,000 K ~ 35,000 K. 

(2) sdOB - Teff is between 35,000 K ~ 40,000 K. 

(3) sdO - Teff is between 40,000 K ~ 80,000 K. 

All of them have helium-burning-cores, and their logarithmic 
surface gravities (log(§)) are between 4.5 and 6.5 (cgs). We dis- 
tinguish sdB and sdOB stars by effective temperature for discus- 
sion, although their differences are also associated with surface 
helium abundance. In fact, sdOB and sdB stars are mixed in the 
region 34,000 - 38,000 K in the T^ff-logig) diagram. 
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Fig. 1. Logarithmic surface gravities versus effective temperature of observed hot subdwarfs (see typical error-bars in the upper left 
side). Open circles and solid triangles are for sdB and sdOB stars observed by Edelmann et al. (120031 ). Open triangles are for sdB 
stars observed by Lisker et al. (120051) . Crosses are for sdO stars observed by Stroeer et al. (2007 ). Dashed line, dot-dashed line and 
dotted line indicate ZAEHB (the Zero- Age Extreme Horizontal Branch, assuming the mass of hot subdwarf is 0.5 M©.), TAEHB 
(the Terminal-Age Extreme Horizontal Branch, assuming the mass of hot subdwarf is 0.5 Mq.) (Dorman et al. 1 19931 ) and He-ZAMS 
(the Helium Zero- Age Main Sequence) (Paczvriski |1971| ). 



3. 1 . The initial parameter space 

We are interested in the ranges of the initial parameters of the 
progenitors of hot subdwarfs through the stable RLOF channel 
of intermediate-mass binaries, i.e. the primary mass (MO, the 
mass ratio {qi), the orbital periods (logPi). Figs. 2 and 3 show 
the initial parameter space of the progenitors of hot subdwarfs in 
conservative case and nonconservative case. 

From Figs. 2 and 3, we see that the number of binary sys- 
tems experiencing stable mass transfer decreases with increasing 
mass ratio qi. When qi >3.0, a large number of binary systems in 
the conservative case have mass transfer rates higher than 10"^ 
Moyr"' (the light grey region). The number of binaries with sta- 
ble mass transfer in the nonconservative case is more than that 
in the conservative case due to mass loss from binaries. The pri- 
mary lying in the light grey region in Figs. 2 and 3 would become 
a hot subdwarf in ~10^yrs after it experiences a stable mass 
transfer with a rate higher than 10"^ Moyr"' which last ~10^ to 
~10^ yrs. Given qi >4.2 in the conservative case (^i >4.5 in the 
nonconservative case), the mass transfer of all systems in our 
models is dynamically unstable, and the stellar evolution code 
breaks down. 

The formation of different types of hot subdwarfs depends 
greatly on the initial parameters (i.e. initial mass, initial orbital 
periods) of binaries, although this relationship becomes weaker 
if qi >4.0. We depict these relationships as follows: 

(i) We can understand the relationship between the formation 
of different types of hot subdwarfs and the initial mass of their 
progenitors. A star with a larger initial mass could leave a larger 
helium core. If the helium core is ignited, the star becomes a hot 
subdwarf with an earlier spectral type. 



(ii) The relationship between the formation of different types 
of hot subdwarfs and the initial orbital periods is more compli- 
cated than the relationship mentioned above. 

From Figs. 2 and 3, for a given initial mass, such as 
Mi =5.O1M0, we can see that, if qi <3.0, it may produce sdB 
stars (sdOB or sdO stars) when the initial orbital periods are 
shorter than 2.1 days (longer than 2.1 days); if qi >3.0, only 
sdO stars will be produced. This is because an MS star fills its 
Roche lobe and this results in the onset of the mass transfer at 
the MS stage. In this case, a star loses more mass because of a 
longer mass transfer time, so that a smaller helium core and a 
thinner envelope would be left, when the star enters the extreme 
horizontal branch (EHB) stage. 

In addition, if qi <1.5, for a given initial mass, such as 
Mi =5.01Mo, we can see that sdB or sdOB stars (sdO stars) 
may also be produced when the initial orbital period is longer 
than 31.6 days (shorter than 31.6 days but longer than 2. 1 days), 
where the mass transfer begins at the FGB stage. In this case, the 
timescale of the mass transfer is very short, as the primary con- 
tracts soon. This leads to less mass lost from the primary, leaving 
an sdB star with a larger helium core and thicker envelope. 

3.2. Tlie Teff-log(g) diagram 

The evolution tracks of some hot subdwarfs with different initial 
parameters are shown in Figure 4. Panels (a) and (b) are for the 
conservative case while Panels (c) and (d) for the nonconserva- 
tive case. 

As seen from Fig. 4, the effective temperatures Teff of hot 
subdwarfs increase with increasing initial mass while their log- 
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Fig. 2. Initial parameters of progenitors of hot subdwarfs in the conservative case. Pluses (black) are for progenitors of sdB stars, 
circles (dark grey) for progenitors of sdOB stars, and crosses (light grey) for progenitors of sdO stars. Panels (1) ~ (6) indicate the 
mass ratio qi of 1.1, 1.5, 2.0, 3.0, 4.0, 4.2, respectively. 
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Fig. 3. Initial parameters of progenitors of hot subdwarfs in the nonconservative case. Pluses (black) are for progenitors of sdB stars, 
circles (dark grey) for progenitors of sdOB stars, and crosses (light grey) for progenitors of sdO stars. Panels (1) ~ (6) indicate the 
mass ratio qi of 1.1, 1.5, 2.0, 3.0, 4.0, 4.5, respectively. 



arithmic surface gravities log(g) decrease. This is because a star 
with higher initial mass could more easily leave a larger helium 
core. For a given initial mass, log(g) of hot subdwarfs decreases 
with increasing initial orbital periods. 

We find that the relationship between Teg and Pi is some- 
what strange both in the conservative and nonconservative case, 
especially for given qi <1.5 and Mi >4.0Mo. As mentioned in 
Section 3.1, this phenomenon is associated with the mass of the 
envelope of hot subdwarfs. In order to display the relationship 
between Teff and Pi clearly, we plot in Fig. 5 the evolution tracks 
of hot subdwarfs with different initial mass and orbital periods 
in the reff-log(^) diagram when qi =1.5. For a given initial mass. 



such as 5.01 Mq, T^s would increase and then it would decline 
with increasing Pi. 

From Fig. 4, we can see that there is a lack of observed 
hot subdwarfs in the upper left region of the Tgff-log(g) dia- 
gram. This is consistent with our models (Note that the evolu- 
tion curves in Fig. 4 are at the boundary for producing hot sub- 
dwarfs.). In addition, if the mass of the helium core of a star is 
too small to be ignited, it will not experience the core-helium- 
burning stage. So, it is difficult to observe hot subdwarfs in the 
lower right region of the 7'eff-log(^) diagram. 

The evolutionary tracks of hot subdwarfs in our models can 
cover positions of observed hot subdwarfs except for a few 
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observed hot subdwarfs with effective temperatures between 
39,000K and 50,000K and logarithmic surface gravities between 
6.1 and 6.4 (cgs). The origin of these hot subdwarfs is not clear. 



3.3. The mass, the envelope mass and the orbital periods of 
hot subdwarfs 

We have obtained the mass and the envelope mass of hot sub- 
dwarfs which depend strongly on the initial parameters of their 
progenitors, i.e. the primary mass (Mu), the orbital periods (Pi) 
and the mass ratio (^i). The mass and the envelope mass of hot 
subdwarfs for the conservative case are shown in Figs. 6 and 7, 
and they are shown in Figs. 8 and 9 for the nonconservative case. 

The relationship between the envelope mass of hot subd- 
warfs and the initial parameters of their progenitors can help 
us to understand the post-evolution of hot subdwarfs. As shown 
by Figs. 7 and 9, the envelope mass has a weak dependence on 
the initial mass ratio if qi >2.0 while the envelope mass has a 
clearly wider range if <1.5. For a given initial mass and ini- 
tial mass ratio, the longer the initial orbital periods, the higher 
the envelope mass. Hot subdwarfs would enter asymptotic giant 
branch manque (AGB-manque) stage after their core helium ex- 
haustion if their envelope masses are between ~0.01 and 0.04 
Mq; hot subdwarfs with an envelope mass of ~ 0.04 ~ 0.08 Mq 
would undergo the early asymptotic giant branch (early-AGB) 
stage; hot subdwarfs with the envelope mass larger than ~0.08 
Mq would become normal asymptotic giant branch (AGB) stars 
(see Fig. 10 for an example). 

As seen from Figs. 6-9, a few sdB stars have progenitors 
with a mass of 3.98 - 5.01 Mq if q^ <3.0, while the majority of 
sdB stars come from stars with mass of 2.23 - 3.55 Mq. This 
is because a primary with a high initial mass would lose more 
mass after a long mass transfer phase due to the primary filling 
its Roche lobe at the MS stage. The progenitors of sdO (sdOB) 
stars have high mass in general, e.g. >5.01 Mq (3.98 Mq ~ 4.47 

Mq). 

In particular, we find that if the mass ratio q-, <1.5, there are 
some sdB or sdOB stars with anomalously high mass and enve- 
lope mass ('anomalous' sdB or sdOB stars, see Table 2.) formed 
though the stable RLOF channel, whose progenitors have a high 
initial mass and a long initial orbital period, such as Mh - 5.01 
Mq, Pi - 40.6 days (see Table 3.). These results are consistent 
with the initial-final mass relation discussed by Han et al. (120001) . 
Chen etal. ( [2002112003] ). 

In order to obtain the mass distribution of hot subdwarfs and 
the final orbital period distribution of the binaries in this work, 
we have plotted Fig. 11 and 12 by executing a simple Monte 
Carlo simulation as mentioned in section 2.2. It is seen that the 
mass of hot subdwarfs are in a wide range from 0.33 to 1.22 Mq 
with a peak near 0.45 Mq. The orbital periods of the binaries are 
also in a wide range from 5.0 to 900.0 days, and peak at around 
120.0 days in the conservative case while the peak at around 30.0 
days in the nonconservative case. 

In Fig. 12, we have also plotted the orbital periods of ob- 
served binary systems containing hot subdwarfs (Morales-Rueda 
et al. 120031 Edelmann et al. |2^ Napiwotzki et al. l2004i) as 
short bars with different colors. A few observed binaries with 
long orbital periods could be formed through the stable RLOF 
channel while the majority of observed binaries with short or- 
bital periods probably would come from the CE ejection channel 
(Han etaL |2U02l ). 



Table 1. The mass and the envelope mass of normal hot subd- 
warf stars. 



type 


Mf 








(Mo) 


conservative case: 






sdB 


0.33-0.67 


0.01-0.05 


sdOB 


0.61-0.86 


0.05-0.10 


sdO 


0.75-1.44 


0.07-0.22 


nonconservative case: 






sdB 


0.32-0.61 


0.01-0.05 


sdOB 


0.60-0.85 


0.05-0.09 


sdO 


0.81-1.36 


0.08-0.19 



Mf = the mass of hot subdwarfs; Me„v = the envelope mass of hot 
subdwarfs. 

Table 2. The mass and the envelope mass of anomalous hot sub- 
dwarfs for mass ratios qi =1.1, 1.5. 



type 


Mf 


-^env 




(Mo) 


(Mo) 


conservative case: 






sdB 


0.72-1.06 


0.07-0.16 


sdOB 


0.98-1.23 


0.11-0.19 


nonconservative case: 






sdB 


0.83-1.13 


0.10-0.16 


sdOB 


0.97-1.17 


0.11-0.16 



Mf = the mass of hot subdwarfs; Me„v = the envelope mass of hot 
subdwarfs. 



3.4. The onset of mass transfer at different stage 

The properties of hot subdwarfs are significantly affected by the 
stage where mass transfer began in their progenitors. When the 
onset of mass transfer is at the main sequence, the primary would 
lose an amount of mass. If its helium core is ignited, the primary 
would become a hot subdwarf with a low mass and a low enve- 
lope mass. Table 3 lists the parameters of hot subdwarfs and their 
progenitors where mass transfer begins at the main sequence, 
Hertzsprung gap or first giant branch in the conservative case 
and the nonconservative case. 

From our calculations, we found that the properties of the 
hot subdwarfs strongly depend on the initial orbital period when 
qi <1.5, but the dependence becomes weaker for hot subdwarfs 
with high mass ratios. Figs. 10 and 13 show the Hertzsprung- 
Russell (HR) diagrams and the Tsff-log{g) diagrams of hot sub- 
dwarfs for a given initial primary mass with a different initial 
mass ratio {qi -l.l and 4.0) when mass transfer begins at differ- 
ent stages in the conservative and nonconservative cases. From 
Figs. 10 and 13, we can see that the onset of mass transfer at 
different stages affects not only the surface gravity of hot sub- 
dwarfs but also their effective temperature given qi =1.1. The 
surface gravity decreases with the increase of the initial orbital 
period. For qi =4.0, the onset of mass transfer at different stages 
almost does not affect the effective temperature of hot subdwarfs, 
but the surface gravity decreases slightly with the initial orbital 
period. 

In the progenitors of hot subdwarfs, the possibility of the on- 
set of mass transfer at the main sequence decreases with mass ra- 
tio qi. In order to understand the mass transfer at different stages 
of the primary with different initial mass, we have plotted inte- 
gral evolutionary tracks of the primary with different mass ratios 
in the HR diagram and mass transfer rate in Fig. 14. From Fig. 
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14, we see that the earlier the mass transfer begins, the lower 
the mass transfer rate. However, mass transfer beginning at an 
early stage tends to leave an sdB star. With an increase of ini- 
tial orbital period, mass transfer rates increase steadily, even up 
to 10""* Moyr"', but the timescale of mass transfer of 10^ yr is 
quite short. Some binaries with high initial mass and long initial 
orbital period would be on the brink of experiencing a delayed 
dynamical instability. 



3.5. Birth rate and number of hot subdwarfs from the stable 
RLOF channel 



Table 4 lists the birth rates of hot subdwarfs formed in various 
formation channels. In the table, the first column denotes the 
metallicity (Z=0.02 for Population I); the second column gives 
^crit, the critical mass ratio for the first stable RLOF on the FGB; 
Columns 3-8 list Galactic birth rates for hot subdwarfs (in 10"^ 
yr"') from the stable RLOF channel of intermediate-mass bina- 
ries (our model), the first stable RLOF channel, the first CE ejec- 
tion channel, the second RLOF channel, the second CE ejection 
channel and the helium WD merger channel. 
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Fig. 8. The mass of hot subdwarfs (Mf) versus the mass of their pogenitors (Mu) in the nonconservative case. Pluses (black) are for 
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Fig. 9. The envelope mass of hot subdwarfs (Mem) versus the mass of their pogenitors (Mu) in the nonconservative case. Pluses 
(black) are for sdB stars, circles (dark grey) for sdOB stars, and crosses (light grey) for sdO stars. Panels (1) ~ (6) indicate the mass 
ratio qi of 1.1, 1.5, 2.0, 3.0, 4.0, 4.5, respectively. 



As Table 4 shows, the predicted birth rate of Population I 
hot subdwarfs from the stable RLOF channel of intermediate- 
mass binaries is 0.004 yr"' for the Galaxy. By taking an effec- 
tive Galactic volume of 5 x 10" pc^ (Zombeck l 19901 1. we give 
an average birth rate per pc^ of 8 x 10"'^ pc"^ yr"'. When con- 
volved with the lifetime of hot subdwarf phase, this rate implies 
a total number of hot subdwarfs from the stable RLOF chan- 
nel of intermediate-mass binaries in the Galaxy of 0.8 x 10^. 
Comparing this result with the result of Han et al. (120031 ), we 
conclude that the stable RLOF channel of intermediate-mass bi- 
naries is important. 



In order to compare our results with observed hot subdwarfs 
in the Ts,ff-\og(g) diagram, we plot a grey diagram as Fig. 15. 
Comparing this figure with Fig. 4, we suppose that some hot sub- 
dwarfs with Teff > 30,000 K and log(g)> 6.0 might be evolved 
sdB stars. 



3.6. The effect of the secondary 

Most hot subdwarfs contain a main sequence (MS) or a subgiant 
star in our models. We discuss the effect of the secondary on 
observations for hot subdwarfs with a companion MS star in this 
section. 
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Fig. 10. The HR diagrams and the Teff-log(g) diagrams of hot subdwarfs where the initial mass (Mu) is 5.01 Mq and the mass ratio 
qiis 1.1. Solid curves are for initial orbital periods of 1 .6 days (sdB stars with envelope mass of 0.04 Mq in left panel, envelope mass 
of 0.03 Mq in right panel), dashed for 2.5 days (sdOB stars with envelope mass of 0.05 Mq), dot-dashed for 15.7 days (sdO stars 
with envelope mass of 0.09 Mq), dash-dot-dot-dot for 31.3 days (left panel, sdOB star with envelope mass of 0.12 Mq). Panels (a) 
and (b) are for the conservative case. Panels (c) and (d) are for the nonconservative case. Asterisks indicate that helium is exhausted 
in its center. 




Fig. 11. The distributions of the masses of hot subdwarfs from the stable RLOF channel in this paper Mf is the mass of hot 
subdwarfs. The solid and dotted line indicate the conservative and nonconservative cases, respectively. 



The luminosity of the secondary (Z =0.02) is fitted by Tout 
et al. ( fT996] l as: 

~ y + M^ + 6M^ + sM^ + ^M^ + i]M'^-^ ' 

where M is the mass of a main sequence star, and a,/3, y, 6, s, rj 
are constants (see Tout et al. 1996). Our calculation indicates 
that if qi <2.0, the luminosity of almost all hot subdwarfs is 
lower than that of their companions. This implies that these hot 
subdwarfs would not be detected by ground-based telescopes, 
although they could be found by space telescopes. If qi >3.0, the 
luminosity of most hot subdwarfs will be larger than that of their 
companions, and this indicates that these hot subdwarfs might be 
observed by ground-based telescope. So, we suggest that a small 



fraction of observed hot subdwarfs could have progenitors with 
qi larger than ~3.0. 

4. Discussion 

Han et al. ( 120021 ) have proposed a binary model for the forma- 
tion of hot subdwarfs, including three channels: the CE ejection 
channel, the stable RLOF channel and the double helium WD 
merger channel. In the CE ejection channel, a primary experi- 
ences dynamically unstable mass transfer at the tip of the FGB 
and leads to the formation of a CE (Paczyhski 11 9761 ) . During 
the spiral-in process in the CE, the orbit of the binary shrinks 
due to the friction between binary and envelope. The CE will 
be ejected if the released orbital energy exceeds the binding en- 
ergy in this process, and this process would leave a very tight 
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Fig. 12. Orbital-period distributions for stable RLOF hot subdwarfs. The solid and dotted line indicate the conservative and noncon- 
servative cases, respectively. The red, blue and green bar along the period axis indicate the orbital periods of observed binaries with 
a companion of unknown type, main sequence and compact object, respectively (Morales-Rueda et al. 120031 Edelmann et al. 120051 
Napiwotzki et al. 2004) . 
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Fig. 13. The HR diagrams and the Teff-log^ diagrams of hot subdwarfs with progenitor's initial mass Mu of 5.01 Mq and initial 
mass ratio qi is 4.0. All solid curves are for initial orbital periods of 5. 1 days (sdO stars) and all dot-dashed for 32.3 days (sdO stars). 
Panels (a) and (b) are for the conservative case. Panels (c) and (d) are for the non-conservative case. Asterisks indicate that helium 
is exhausted in its center. 



binary containing the degenerate core of the giant. If the degen- 
erate core experiences a helium flash, the remnant core of the 
giant may be ignited (Castellani & Castellani I1993I I and hence 
appear as a core-helium-burning hot subdwarf in binaries with 
short orbital periods in the range of 0.05 - 40 days (Han et al. 
I2OO2I 1. In stable RLOF channel, a primary undergoes a stable 
mass transfer at the FGB stage. If its helium core was ignited 
after the RLOF, the primary would become a hot subdwarf in bi- 
naries with long orbital periods in the range of 400 - 1500 days 
(Han et al. |2002| l. Our results indicate that the mass transfer can 
also begin at the MS stage (case A mass transfer) and HG stage 
(case B mass transfer). This leads to hot subdwarfs in binaries 
with a wide orbital period range of 5 - 900 days. Compared with 
the CE ejection channel, we found that there are a large number 
of hot subdwarfs in binaries with long orbital periods formed 



through the stable RLOF channel. This would sustain the obser- 
vations by Green, Liebert & Saffer (2000J and Morales-Rueda 
et al. ( I2003I I, who argued that some sdB stars could appear to 
be members of long period binaries. Moreover, it is possible that 
some hot subdwarfs in binaries with a very long orbital period 
could be misunderstood as single stars. 

In general, sdB stars are believed to have a mass of around 
0.5 Mq and a thin hydrogen rich envelope around 0.02 Mq 
(Heber |T986l Saffer et al. lT^ . Han et al. (2002 ) suggested that 
the mass of sdB stars is in the range of 0.33 - 0.68 Mq. In our 
models, the mass of most sdB stars is in the range of 0.32 - 0.67 
Mq with a thin envelope of 0.01 - 0.05 Mq, but a few sdB stars 
have an anomalously high mass of 0.72 - 1.13 Mq, and a thick 
envelope of 0.07 - 0.16 Mq. The reasons have been described in 
section 3. The mass of sdO stars is from 0.75 Mq to 1.44 Mq 
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Fig. 14. Integral evolutionary tracks (top panels) and mass transfer rate (bottom panels) of the primary with initial mass of Mn =5.01 
Mq in the non-conservative case. The left panels are for q -I. I, and the right panels are for q -4.0. In the left panels, solid curves 
are for the initial orbital periods of 1.6 days, mass transfer beginning at the MS stage (sdB star, with envelope mass of 0.03 Mq), 
dashed for 15.7 days, mass transfer beginning at the HG stage (sdO star with envelope mass of 0.09 Mq) and dot-dashed for 43.6 
days, mass transfer beginning at the FGB stage (sdB star with an envelope mass of 0.14 Mq). In the right panels, solid curves are 
for the initial orbital periods of 5.1 days, mass transfer beginning at the HG stage (sdO star with an envelope mass of 0.08 Mq) and 
dashed curves for 32.3 days, mass transfer beginning at the HG stage (sdO star with an envelope mass of 0.09 Mq). 




Fig. 15. Comparison of our results to the atmospheric parameters of sdB stars (open circles Edelmann et al. 12003 1, open triangles 
Lisker et al. 2005), sdOB stars (solid triangles Edelmann et al.'2003|l, sdO stars (crosses Stroeer et al. 200^. The largest cross is for 
HE 1115-0631 (orbital period: 5.87 d, Napiwotzki et al. 2004). Our results are shown as shaded reff-log(g)-boxes, where a higher 
hot subdwarf density per box corresponds to darker shading. The grey scale is shown in the right of the figure. 



and its envelope mass is higher than 0.07 Mq, even up to 0.22 
Mq. The mass and the envelope mass of sdOB stars are between 
sdB stars and sdO stars. In addition, Han et al. (2000) and Chen 
et al. (I2OO2I 1 discussed low- and intermediate-mass close binary 
evolution and the initial-final mass relation of these binaries. By 
calculating a binary evolution sequence with an initial mass ra- 
tio between 1 . 1 and 4.0, and the onset of the stable RLOF at the 
Hertzsprung gap stage, they found that the remnant mass of the 



primary depends on the initial mass ratio or orbital period for a 
given primary mass. These results are similar to ours, but binary 
evolution calculations with the onset of mass transfer at the main 
sequence stage have been included in our models. 

Elson et al. (I1998I I discovered a blue star with Jeff ~ 31,500 
K, log{g)x 4.4, logL/Lo ~ 3.0 (Burleigh et al. [19991 ) in the 
young cluster NGC 1818 in the Large Magellanic Cloud (LMC), 
which has a main-sequence turnoff mass of ~7.5 - 9.5 Mq, an age 
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Table 3. Parameters of some typical hot subdwarfs and their pro- 
genitors. 



Mii / Mf(type) 


9i 


P-JPi 


-^env 


'age 


(Mo) 




(days) 


(Mo) 


(yrs) 


conservative case: 










3.55/0.38(sdB) 


1.1 


1.5/151.1 


0.02 


3.8x10** 


3.55/0.56(sdB) 


1.1 


3.5 / 116.1 


0.04 


0.9x10* 


3.55/0.53(sdB) 


4.0 


3.5/11.5 


0.04 


1.1x10* 


5.01/0.54(sdB) 


1.1 


1.6/ 147.8 


0.04 


1.0x10* 


5.01/0.65(sdOB) 


1.1 


2.5 / 137.3 


0.05 


0.6x10* 


5.01/0.93(sdO) 


1.1 


15.7/317.2 


0.09 


0.2x10* 


5.01 / l.Ol(sdOB) 


1.1 


31.3/483.8 


0.12 


0.2x10* 


5.01 / 1.06(sdB) 


1.1 


40.6/545.5 


0.16 


0.1x10* 


5.01 /0.84(sdO) 


4.0 


5.1/13.1 


0.08 


0.2x10* 


5.01 /0.91(sdO) 


4.0 


32.3/60.8 


0.09 


0.2x10* 


non-conservative case: 










3.55/0.33(sdB) 


1.1 


1.5/27.9 


0.01 


7.9x10* 


3.55 / 0.54(sdB) 


1.1 


3.5 / 29.4 


0.04 


1.0x10 


3.55/0.53(sdB) 


4.0 


3.5/6.9 


0.04 


1.1x10* 


5.01 /0.46(sdB) 


1.1 


1.6/27.2 


0.03 


1.7x10* 


5.01 /0.61(sdOB) 


1.1 


2.5/29.0 


0.05 


0.7x10* 


5.01/0.91(sdO) 


1.1 


15.7 / 100.5 


0.09 


0.2x10* 


5.01/0.98(sdOB) 


1.1 


41.6/483.2 


0.12 


0.2x10* 


5.01 / l.OO(sdB) 


1.1 


43.6/526.0 


0.14 


0.1x10* 


5.01 / 0.84(sdO) 


4.0 


5.1/9.1 


0.08 


0.3x10* 


5.01/0.91(sdO) 


4.0 


32.3/59.4 


0.09 


0.2x10* 



M\i = initial mass; Mf = the mass of hot subdwarfs; q\ = the mass 
ratio; P, = initial orbital period; Pf = orbital period of hot subdwarf 
binaries; Me„v = the envelope mass of hot subdwarfs; Z^gc = the lifetime 
of a hot subdwarf. More parameters can be seen from Figs. 2, 3, and 6 - 
9. 



Table 4. Birth rates of hot subdwarfs from the stable RLOF 
channel of intermediate-mass binaries (our model) and other 
channels (Han et al. l2003l l (in 10"^ yr '). 



z 




ViMR 


vfr 


vfce 




VSCE 




0.02 


1.5 


4.18 


29.89 


8.41 


0.00 


8.38 


17.22 


0.02 


1.2 


4.14 


22.25 


10.71 


0.00 


5.43 


7.99 



Z=metallicity; g„„=the critical mass ratio for the first stable RLOF 
on the FGB; v;A},;=Galactic birth rates of hot subdwarfs from the sta- 
ble RLOF channel of intermediate-mass binaries; VfR=the first stable 
RLOF channel; VfC£=the first CE ejection channel; yi7;=the second 
RLOF channel; vscE=tiie second CE ejection channel; yM=the helium 
WD merger channel. 



of ~2 - 4 X 10^ yr ( Will et al. [T995] l and a metallicity of -0.005 
(Kerber & Santiago |2005]l. This star was considered as a candi- 
date of luminous WD. However, Liebert (119991 1 and Burleigh et 
al. (11999 ) pointed out that the star would be a post-EHB star or 
EHB star, rather than a WD, and it would lie in the Galactic halo 
or the disc of the LMC instead of cluster NGC 1818. If this blue 
star was in a binary system, we suggest that the star could be in 
the cluster NGC 1818 and form through the stable RLOF channel 
of intermediate-mass binaries, but this needs more observational 
and theoretical evidence to confirm. We consider that hot sub- 
dwarfs could exist in young clusters, and some of them might 
have high mass (> 0.7 Mq). We will focus on hot subdwarfs in 



young clusters and calculate a binary evolution sequence with 
different metallicity (Z) to study their properties, formation and 
evolution in the future. 

Mass transfer efficiency is an important factor aff'ecting the 
initial parameter spaces, as lower mass transfer efficiency will 
lead to more binaries experiencing stable RLOF rather than the 
formation of a common envelope during rapid mass transfer. Past 
work (Paczytiski & Ziolkowski 1967^ Refsdal, Roth & Weigert 
[T9741 De Greve & De Loore [19921 Chen & Han |2002l [2^31 1 
employed mass transfer efficiency j3 - 0.5 for the calculation of 
stellar evolution in nonconservative cases. However, De Mink et 
al. ( 120071 suggested that mass transfer efficiency should not be 
a single constant by comparing their models with a sample of 50 
double-lined eclipsing binaries in the Small Magellanic Cloud, 
which could be affected by spin up of the accreting star and tidal 
interaction. They also found that initially wider systems tend to 
favor less conservative models, since accreting angular momen- 
tum will speed up the rotation of the accreting star, resulting 
in mass loss along its equator region, as discussed by Wellstein 
(1200 11 and Petrovic et al. ( 120051 1. Our results are in line with this 
conclusion. Due to the uncertain of the mass transfer efficiency 
(De Mink et al. 120071 1. the birth rate and the total number from 
our models need further data to constrain. However, it is reason- 
able to estimate that a high birth rate and large total number of 
hot subdwarfs via stable RLOF channel will suggest less conser- 
vative models. 

Other effects on the birth rate and total number of hot sub- 
dwarfs in our models would arise from the life-time of the sec- 
ondary, which will become too short, such that it is comparable 
to or even smaller than the lifetime of the hot subdwarf phase 
(typical value: 2.0x10** yr) in some extreme cases. In these two 
cases, the secondary presumably evolved rapidly to the giant 
branch or even to later stages after fast mass transfer, while the 
primary could be a giant or a hot subdwarf or a helium star, 
which could be associated with a 'post-Algol' binary (Nelson 
& Eggleton 2001 )■ As a result, we could overestimate the birth 
rate and the total number of hot subdwarfs via the stable RLOF 
channel. 

5. Conclusions 

We present the properties of hot subdwarfs from the stable 
RLOF channel of intermediate-mass binaries by computing pop- 
ulation I binary evolution sequences both in conservative and 
nonconservative cases, where the primary has an initial mass in 
the range of 2.2 - 6.3 Mq, initial mass ratio qi =1.1, 1.5, 2.0, 
3.0, 4.0, 4.2, 4.5, and onset of RLOF at the MS, HG, and FGB 
stage. Due to the effect of the secondary, the ratio of the hot sub- 
dwarf bianries which might be observed by ground-based tele- 
scope would be larger around 3.0. We summarize our results as 
follows: 

1 . The birth rate of hot subdwarfs from the stable RLOF chan- 
nel of intermediate-mass binaries is ~4.2xl0"'' yr"', which 
is similar to the first CE ejection channel. However, the cur- 
rent observations indicate that only a few hot subdwarfs bi- 
naries could be formed via the stable RLOF channel (Fig. 
12). This could be associated with the rotation of the accret- 
ing star and the evolution of the secondary. 

2. We have obtained the initial parameter space of the progeni- 
tor binaries of hot subdwarfs from the stable RLOF channel 
of intermediate-mass binaries (Figs. 2 and 3). 

3. Our results indicate that the hot subdwarfs from the stable 
RLOF channel of intermediate-mass binaries have a mass in 
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the range of 0.32 - 1.44 with a wide envelope mass of 
0.01 - 0.22 Mq. Most sdB stars have a mass in the range of 
0.32 - 0.67 Mq with an envelope mass of 0.01 - 0.05 Mq. 
Most sdOB stars have a mass in the range of 0.60 - 0.86 Mq 
with an envelope mass of 0.05 - 0.10 Mq. SdO stars have a 
mass in the range of 0.75 - 1.44 Mq with an envelope mass 
of 0.07 - 0.22 Mq. Furthermore, we found that some sdB 
or sdOB stars have anomalously high mass, in the range of 
0.72 - 1.13 Mq and 0.97 - 1.23 Mq, respectively. The enve- 
lope mass is in the range of 0.07 - 0.16 Mq and 0.11 - 0.19 
Mq, respectively. We expect that observations could confirm 
these theoretical results in the future. 

4. The orbital periods of the binaries have a wide range of 5.0 
- 900.0 days, while the peak is at around 120 days for the 
conservative case and 30 days for the nonconservative case. 

5. Our results favor that hot subdwarfs, formed through the 
RLOF channel of intermediate-mass binaries, will be found 
in young clusters in the future. 
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